Full counting statistics of strongly non-Ohmic transport through single molecules 
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We study analytically the full counting statistics of charge transport through single molecules, 
strongly coupled to a weakly damped vibrational mode. The specifics of transport in this regime - a 
hierarchical sequence of avalanches of transferred charges, interrupted by "quiet" periods - make the 
counting statistics strongly non-Gaussian. We support our findings for the counting statistics as well 
as for the frequency-dependent noise power by numerical simulations, finding excellent agreement. 
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Introduction. — A prime qualitative difference of trans- 
port through single molecules as compared to artificial 
nanostructures lies in the role of the vibrational motion 
of the nuclei. This aspect is at the focus of current ex- 
periments 0,0,0, and is also beingstudied theoretically 
within a number of approaches 0, H, 0, 0, IE E3 Tne 
incorporation of molecular vibrations (phonons) into the 
theoretical description is mostly done within simplified 
(phcnomenological) models, as opposed to purely elec- 
tronic first-principles studies [TH . 

Even within minimal models involving one molecular 
orbital coupled to a single vibrational mode, "unidirec- 
tional" transport (i.e., for voltages large compared to 
temperature) depends radically on the strength of the 
electron-phonon coupling, already at the qualitative level 
[E IHj- F° r weak and intermediate coupling [l3| . 
transport is adequately described in terms of individual 
electron transitions. By contrast when vibrational relax- 
ation is slow, transport in the regime of strong electron- 
phonon coupling is appropriately captured within a sce- 
nario of avalanches of transferred electrons, with expo- 
nential spreads of height and duration |lfjj . 

More quantitatively, the time dependence of the cur- 
rent in the strong-coupling regime can be presented as 



I(T)=ri°\T-t 1 ) + f£"(T-t 1 -t 



(0), 



(1) 



where U are the time intervals between avalanches (quiet 
periods). These intervals are much longer than the typ- 
ical duration of an avalanche which occurs dur- 
ing the sparse periods when the vibrations are excited. 
The random function (t) (which is nonzero only for 
times |t| < t^ 3 ') describes the passage of a large number 

J drfi (t) = Ni 3> 1 of electrons through the molecule 
during the ith. avalanche. Moreover, a numerical study 
of the avalanches revealed their self-similar hierar- 
chical structure, see Fig. 1. Quantitatively, this struc- 
ture manifests itself in the fact that, during the time of 
an avalanche ~ t^°\ each function itself takes the 
form of Eq. QJ, with replaced by random functions 
fj , which describe avalanches of the first generation in- 



terrupted by quiet periods. Again, these quiet periods 
are much longer than the characteristic time scale rW of 
the functions . For times shorter than rW, the func- 
tions fj have the form of Eq. (JJJ with corresponding 

(2) 

second-generation avalanches, f„ , having even shorter 
time-scale, t^ 2 \ and so on jl4| . Numerical results sup- 
porting this scenario are shown in Fig. 1. 

The above discussion implies that the statistical prop- 
erties of charge transport through a molecule in the 
regime of strong electron-phonon coupling and through a 
conventional nanostructure are drastically different. For 
a nanostructure, all f\ are 8 functions, so that Ni = 1. 
Hence, the distribution function Pt(Q) of the net trans- 
mitted charge Q during time T (full counting statistics 
[lfij ) is completely encoded in the distribution of the wait- 
ing times ti for single-electron transitions. This distribu- 
tion reflects the details of the transport mechanism, and 
might be quite nontrivial ^3- Nevertheless, with all ti 
being of the same order, the full counting statistics differs 
only weakly from a Gaussian distribution. Small devia- 
tions are caused by correlations [Til Hp . interactions [l9j . 
or the influence of the environment |20| , and have been 
extensively studied theoretically. 

By contrast, the counting statistics of avalanche- type 
transport is insensitive to the details of the passage of 
a single electron through the molecule, since the num- 
ber of electrons involved in each avalanche is large. In- 
stead, the counting statistics is governed exclusively by 
the transition rates between different vibrational states. 
These rates have a simple structure in the limit of strong 
coupling which allows us to develop a complete analyt- 
ical theory for the regime of avalanche- type transport. 
In particular, we demonstrate in this paper that the full 
counting statistics Pt{Q) is given by a concise analytical 
expression, which is strongly skewed at "short" times (~ 
zero-order quiet period) and evolves into a Gaussian only 
for very large T. Along with the counting statistics, we 
also study analytically how the hierarchy of avalanches 
manifests itself in the frequency dependence of the noise 
power S(uj). Our analytical results are in excellent agree- 
ment with numerical Monte-Carlo (MC) simulations. 
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FIG. 1: Hierarchical character of transport. Left: Three gen- 
erations of self-similar MC plots for A = 4 and eV = 3hu> v , 
showing the net-transferred charge TV and phonon state q as 
functions of time. (ui v : phonon frequency; p: density of states 
of the leads; to: molecule-lead coupling). Right: Comparison 
of the fixed-point distribution for the transferred charge per 
generation-g avalanche to numerical simulations for q = 0, 1, 2 
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Full counting statistics. — Since different zeroth- 
generation avalanches are statistically independent, it is 
easy to derive a relation between the counting statistics 
Pt(Q) of the net charge Q and the conventional counting 
statistics </?t(^) of the number of zeroth-generation 
avalanches n during the time interval T . Indeed, using 
the definition Pt{Q) = (S(Q - Y!j=i N j))Nj,n, and a 
Fourier representation of the RHS, one obtains 

da A „r> 



MQ) - I ^ 



r ~ i " 

[V (a)\ tpr(n), (2) 



where Vo( a ) = (exp (— iaNj))^. denotes the Fourier 
transform of the distribution function 'Pq(N) of the total 
charge passing per zeroth-generation avalanche. The du- 
rations of the quiet periods obey Poisson statistics so that 
(fiT( n ) = exp(— n^rffp/nl. Here, Ht denotes the average 
number of zeroth-generation avalanches within time T. 
Substituting this form into Eq. © and performing the 
summation over n yields an expression for the counting 
statistics similar to the Holtsmark distribution 12 ll . 



Pt{Q)= I ^expjiaQ 



11T 



P (a) - 1 



(3) 



Thus, the problem of the counting statistics is reduced 
to finding the distribution Vq(N). Two facts allow us 
to find Vq(N), namely (i) the self-similar structure of 
avalanches and (ii) the large number n q of generation- 
(q + 1) avalanches within a given generation-g avalanche. 

Our basic observation is that we can derive a recur- 
sion relation, relating the distribution functions V q (N) 
and , P g +i(TV) of the total passing charge (TV^ and 



N^ q+1 \ respectively) per avalanche for subsequent gen- 
erations. This recursion follows from the obvious facts 
that AT(9) = YTjli TVj 9+1) and that different avalanches 
of a given generation are statistically independent. By 
analogy with the derivation of Eq. © , we thus obtain 



r - 1 71 

\P q+ x{a)\ Pq (n), (4) 



where p q (n) denotes the distribution function of n q . To 
proceed further one has to specify the form of the distri- 
bution p q (n) . This distribution is governed by the micro- 
scopic characteristics of the Franck-Condon transitions. 
We demonstrate below that p q (n) = (l/n q ) exp(— n/n q ). 
Upon substituting this form into Eq. Q), the summa- 
tion over n on the RHS can be easily performed and we 
obtain, after a Fourier transform of both sides, 



1 



V q (a) = - 



n q l-V q+1 (a)exp(-l/n q ) 



(5) 



The distribution V q can now be obtained from this equa- 
tion by writing its general solution as V q (a) — [1 + 

iaN (q) + c q (a~N {q) ) 2 + . . Inserting this into Eq. ©, 
we find that the numerical coefficients c q flow to zero 
with q by virtue of the small parameter l/n q . Thus, 

the solution V q (a) = [1 + iaN^]^ 1 with Fourier trans- 
form V q {N) = 6(N)exp(-N/N {q) )/N (q) can be viewed 
as a fixed point of the recursion equation Eq. (0J and 

since TV 9 = n q N 9+1 , self-similarity is obeyed asymp- 
totically. The existence of this fixed-point solution can be 
viewed as a consequence of remark (j) which implies that 
up to rescalings, the distribution functions V q (N) have 
the same functional form for all q. Fig. ^ numerically 
confirms this result for three different generations. 

With P q (N) established, we obtain the counting statis- 
tics by substituting Vq{o) = (1 + iaN^y 1 into Eq. © 
and performing the integral. This yields 



Pt(Q) 



Q 
v< ri > 




Here I\{z) denotes a modified Bessel function. Eq. © 
is our central result. It is nicely confirmed by our MC 
results, as shown in Fig. [21 and describes the evolution of 
the counting statistics between the following two trans- 
parent limits, (i) Short times, ut = T/(ti) <C 1: Using 
the expansion I\(z) « z/2 for z <C 1 we obtain from Eq. 
© 



Pt(Q) 



"t 6(Q) + (n T /N i0) )e- Q / Nm 



(7) 



Typically only a few electrons are transmitted through a 
molecule. The long tail described by the second term in 
Eq. Q arises from realizations where an avalanche occurs 
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ln(T/T ) = 13.5 
Ht = 1.78 
ln(r/T ) = 14.9 
nr = 7.15 
ln(T/T ) = 15.6 
n T = 14.3 
ln(T/T ) = 17.7 
n T = 114 




FIG. 2: (Color online) Evolution of full counting statistics 
Pt(Q) for four different time intervals T with A = 4.0 and 
eV — 3huj v . The MC data (solid lines) are in excellent 
agreement with the analytical full counting statistics, Eq. (|SJ, 

(dashed lines), with N {0) = 91.2 and n T = 2.4 ■ Kr 6 T/T (no 
fit), as well as To = (2tt pt 2 /h)' 1 . 



within the time T and reflects the spread of charge within 
a single avalanche, (ii) Long times, Ht 3> 1: Upon sub- 
stituting the large- z asymptote of h(z) into Eq. |JB}, it 
is easy to see that the second term has a sharp max- 
imum centered at Q = n^N^, which is the average 
charge passed through the molecule after a large num- 
ber of avalanches. Expansion of the exponent around 
the maximum yields the Gaussian 



P T (Q) ~ (V2 7 r f j Q )- i exp 



-(Q-n T N i0) 



7(0) 



(8) 



1/2" v 

with a width o~q = (2t%t) N . This width is twice 
the width expected from the fluctuations of the waiting 
times. This enhanced broadening is due to fluctuations 
of the charge passed per avalanche. These additional 



fluctuations also manifest themselves in the noise char- 
acteristics of transport as analyzed below. 

Microscopic derivation of p q (n). — The distribution 
p q {n) is obtained by averaging the Poisson distribution of 
n for a given avalanche duration over the distribution of 
durations. On microscopic grounds, the latter distribu- 
tion is a simple exponent, which immediately transforms 
into p q {n) — (1/Hq) exp(—n/n q ) since n q is large. 

To see this, we note that the duration of a 

generation-g avalanche is determined by the waiting 
times in the vibrational state q + 1 since the dura- 
tions of intermittent higher-generation avalanches can 
be neglected. Two processes terminate a generation- g 
avalanche: a direct transition from q + 1 to q or a transi- 
tion back to q during a generation- (g+ 1) avalanche. De- 
noting the total rate for both processes by T q , we obtain 
an exponential distribution of durations T q exp(— T^r^). 

Noise spectrum S(u>) of avalanche-type transport. — We 
first derive a general expression for S(u) assuming arbi- 
trary distributions Vo(N) of the avalanche magnitudes 
and W(t) of the waiting times. For frequencies smaller 
than 1/t(°', we have fP(t) ~ Ni5(t) in Eq. Q. Using 
Fourier representations of the S functions and averaging 
over the U and Ni, the average current becomes 



(I(T)) - (Ni] 



da 
2tt' 



iaT 



W(a) 
1 - W{a) 



(9) 



where W(a) = (exp(—iati))t t denotes the Fourier trans- 
form of W(t). In the long-time limit, only small values 
of a contribute to the integral in Eq. (JSJ so that we can 
use the expansion W(a) = 1 — ia(U) — {l/2)a 2 (tj) + . . .. 
Inserting this expansion into Eq. keeping only the 
leading order in a and performing the contour integration 
over a, we recover the obvious result (I(T)) = (Ni) / (ti). 
Similarly, we express the current-current correlator as 



(KTi)I(T 2 )) ■ 

da dp 
2tt 27r f 



dadP_ e _ ioli 
2tt 2tt 



i-^/3T 2 Aj Vie iat 1 + N2e Mt!+t 2 ) + _ _ ][ Nie Wi + N2e iP(H+t2) +...]\ 



r!-i0T 2 ) (Nj) 2 W(a + f3) 



1 



l-W{a + /3) \l~W(a) 1-W(J3) 



-1 + 



{{N f)-(N^)W(a + (3) 
1 - W(a + 13) 



(10) 



The last equality follows upon tcrm-by-term averaging 
and resummation of the series. To access the limit of long 
times T = (Tx + T 2 )/2, we introduce w = (a - /3)/2 and 
Q = a + p. Then, the exponent in the integrand in Eq. 
(|10(l assumes the form exp(iu)T — iflT) with r = T% — T±. 
The limit T — > oo can now be taken in analogy with the 
derivation of (I(T)) above. The integrand can be directly 



identified with the noise spectrum S(u>), so that 



SM = ^ (N^ 



1 - W(u) 1 - W{-lu) 
+ ((Nf) (N^) 



(11) 



Taking the zero-frequency limit requires one to keep 
terms of order uj 2 in the expansion of W(±ui). In this 
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way, the Fano factor F = S(lu — 0)/2e(I) becomes 



F = (Ni) 



(if)-(ti) 2 , (N?)-(Niy< 



(ti 



(Ni, 



(12) 



This equation allows for a transparent interpretation: 
Noise originates from two sources, namely the fluctua- 
tions in the intervals between avalanches and the fluctu- 
ations in the transmitted charge per avalanche. In the 
conventional situation where Ni — 1 for all i, the Fano 
factor is given by the fluctuations of the waiting times 
ti for a transition in which an electron passes either di- 
rectly or sequentially from the left to the right lead. For 
example, for transport through a symmetric junction in 
the Coulomb-blockade regime, one immediately recovers 
F = 5/9 [22J when taking into account that the rates 
of entering and leaving the dot are related as 2:1 due to 
spin. For the specific distributions adopted in our model, 
both terms in Eq. i|12[l contribute equally, and the Fano 

factor reduces to F — 2N , which, in agreement with 
Eq. (jSJ , is twice the value expected for a fixed magnitude 
of avalanches. This is confirmed by numerical results. 

For frequencies larger than 1/t^°\ the "fine structure" 
of the avalanches described by the functions in Eq. 
l|T)l must be taken into account. This fine structure 
can be incorporated into the noise spectrum Eq. 1111) 
by replacing (A^) 2 by (/>))</(/3)> and (N?) - (N) 2 
by (/>)/>)} - (/(«)>(/(/?)>, where /(a) denotes the 
Fourier transform. Explicitly employing the Poisson dis- 
tribution of the waiting times leads to the remarkable 
simplification [1 - IF(cj)]- 1 + [1 - W(-u)}- 1 -1 = 1. In 
this way, we obtain 



S(w) 



-</(w)/(-w)>. 



(13) 



For frequencies of order u> ~ 1/tq (where r q denotes the 
average waiting time (ti) at level q of the hierarchy), we 
can ignore the fine structure of the avalanche and replace 



ivj 0) 



Thus, we find S(u) = 2([N ( i 0) } 2 ) /t . At 



higher frequencies u) ~ 1/ti, the function f(t) is resolved 
into avalanches of generation q = 1 . Then, we can write 
= / dT J dre*"(f(T + r/2)/(T - r/2)). 
Up to the integral over T, this expression is analogous to 
S(lo) itself, with zeroth-gencration quantities replaced by 
corresponding q = 1 quantities. For the frequencies of in- 
terest, we therefore find S(w) = (2/t )(t(°) ([N^] 2 ) / Tl ). 



T\riQ and N 



Using the obvious relations 
rioN^ ' and generalizing to arbitrary q, we find 



(0) 



Sq+l — 



Sn . 



(14) 



Here, we define S q — S(uj ~ l/r g ) so that Eq. ifHjl pro- 
vides a rule for extending the noise spectrum to progres- 
sively higher frequencies. 



The essential microscopic inputs are the ratios T q+ i/r q 

and N^ q+1 ^ / N^ . Both ratios are determined by overlaps 
of displaced vibrational wavefunctions 10] . The rate 1 /r q 
is dominated by the transition q —* q+l. Thus, it involves 
the overlap of neighboring harmonic oscillator states. By 
contrast, is inversely proportional to the transition 
rate from a highly excited phonon levels to the gth vi- 
brational level. The difference between these two rates is 
thus that the first involves four wavefunctions with index 
of order q, while the second involves only two. As a re- 
sult, we can immediately establish from a quasiclassical 

evaluation of the matrix elements that T q /(N < ' q ' > ) 2 is es- 
sentially independent of q. With this input, we conclude 
that S(uj) ~ uj~ a with exponent a = 1/2. Since the noise 
power does not depend sensitively on oj in finite intervals 
around 1/t 9 , this power law should be superimposed with 
steplike features in S(u>). These conclusions agree with 
numerical simulations (see Ref. 01) over several orders 
of magnitude in frequency. 

Conclusions. — Our complete analytical description for 
the full-counting statistics and the frequency-dependent 
noise power of self-similar avalanche-type transport was 
made possible by the fact that current flow is essen- 
tially unidirectional. We emphasize that our arguments 
are quite general, with rather limited microscopic in- 
put, making our results potentially applicable far be- 
yond the particular realization [10| of avalanche-type 
transport considered in the numerical simulations. Fi- 
nally, we remark that direct vibrational relaxation (with 
rate 7 rc i), neglected so far, only gradually suppresses 
avalanche- type transport. Indeed, one readily argues 
that ~ l/7rei over a wide range of relaxation times, 
leading to S(0) ~ l/ 7r 2 ol . 
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